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1 Dirty-Fermi-liquids 



It is generally expected that in the absence of strong external fields and 
when lattice effects can be ignored, weakly or even moderately interacting 
systems will remain a Fermi- liquid at sufficiently low temperatures. This 
simple expectation, we now understand, is wrong in two dimensions as even 
the weakest disorder has been shown to destabilize the Fermi-liquid at very 
low temperatures [1, 2], driving the system into a strongly-correlated unusual 
metal or into an insulator [3] depending on the initial (high temperature) 
resistance of the system. Indications are that the the two phases are kept 
apart by a quantum critical point (QCP) at a finite critical resistance [3] . 
We argue that the observed metal-insulator transition (MIT) [4] in two- 
dimensional electronic systems is a manifestation of this break down [5, 6]. 

We are unfortunately very far from understanding the nature of the ul- 
timate metallic state at absolute zero. We do, however, understand the pro- 
cess by which the Fermi- liquid is destabilized. In the diffusion regime, i.e., 
for temperatures (frequencies) less than the elastic scattering rate, T < l/r, 
diffusion modes and not quasi-particles are the low lying propagating modes. 
Observable effects of quantum diffusion accumulate over long distances and 
long times and thereby develop only at much lower temperatures T < \e~" , 
where cr is the dimensionless conductance of the sample. Therefore, non- 
trivial effects are most clearly visible only in the vicinity of the transition 
where cr ~ 1. To follow the effects over a wide range of temperature requires 
a reasonable value for l/r, putting practical constraints on the observability 
of quantum diffusion effects. The scattering rates are much too low in ultra- 
high mobility modulation doped or HIGFET structures (l/r <^ 500mK); 
silicon MOSFETs, on the other hand, have mobilities that are not too high 
leading to l/r 1 - 2 K in the vicinity of the MIT, making them ideal 
playgrounds for studying quantum diffusion. It should be noted that in the 
absence of any other correlation-induced phase transition, all 2D systems 
will at sufficiently low temperatures be destabilized. 

The theory of the dirty-Fermi-liquid which details the destruction of 
the Fermi-liquid phase up to and including the quantum critical point is 
reviewed here. 

1.1 Experimental systems 

The mobility of 2D quantum wells vary over orders of magnitude, rang- 
ing from ~ 20, 000 - 40, 000 cmVVs as in typical n-type Si-MOSFETs, to 
hundreds of thousands and even millions as in ultra-high mobility p-type 
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GaAs heterostructures. The source and therefore the nature of the impurity 
potentials vary greatly between these systems. It is, therefore, important 
to separate the effects of long versus short ranged random potential fluc- 
tuations on the scale of the interparticle distance, especially at very low 
carrier concentrations when the wavelength of the carriers is considerably 
increased. Screening of the impurity potentials becomes difficult in these 
ultra-high mobility systems; the resulting non-linear screening can induce 
non-trivial temperature dependencies in the transport properties of 2D sys- 
tems and can in the extreme situation also lead to domain formations of 
high and low density regions [7] preventing us from experimentally studying 
the fundamental problem of the effects of random quantum scattering in a 
Fermi-liquid, whose full solution is unknown. 

We argue below that MOSFETs provide an ideal environment to study 
the properties of a 2D dirty- Fermi- liquid. The moderate values of the mo- 
bility of these systems results in an impurity elastic scattering rate slightly 
larger than a Kelvin. Comparing with the Drude conductance a = n/i, 
expressed in units of e^//i, where n is the electron density and /i the mo- 
bility, moderately high mobility implies that the density in the vicinity of 
the quantum conductance is of the order of 10^^ cm~^. The correspond- 
ing Fermi-temperature is of the order of 10 K, and is therefore degenerate 
already at not so low temperatures. Also, because the electron-phonon cou- 
pling in silicon is weak due to the presence of inversion symmetry in the bulk, 
electron-phonon scattering contribution to the resistivity is negligible at low 
temperatures. Finally, another important characteristic of MOSFETs, in 
strong contrast to heterostructures, is that the elastic single particle scat- 
tering rate is comparable to the transport scattering rate; they typically 
differ by a factor of 10 in other systems. Quantum scattering therefore is 
the dominant scattering mechanism in MOSFETs. 

The Coulomb interaction energy at the densities of interest are typically 
10 times larger than the Fermi-energy, i.e., = Ec/Ep ~ 10, where Ec is 
the Coulomb energy of neighboring electrons. At these r<j values, detailed 
numerics indicate no phase transitions, at least up to ~ 25 where a 
paramagent to ferromagnet transition is found to occur in a system with 
single valley (no transition has been found in a two- valley system) [8, 9]. At 
still higher ~ 35, the electron liquid is expected to form a Wigner crystal, 
although none has been observed experimentally [10]. These studies suggest 
a Fermi-liquid with moderately strongly interacting quasi-particle scattering 
amplitudes around 10. 

MOSFETs have the additional advantage that they are a multi-valley 
system. The conduction band of an n-(OOl) silicon inversion layer has two 
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almost degenerate valleys located close to the X-points in the Brillouin 
zone [11]. The presence of two equivalent valleys further enhances the cor- 
relation effects [5, 6] making MOSFETs ideal systems to study the effects of 
quantum disorder scattering in a Fermi-liquid. In the following we consider 
the number of equivalent valleys to be equal to A^. 

2 Diffusion regime 

Fermi-liquid renormalizations, like mass, g-factor, or equivalently, specific 
heat and spin-susceptibility, occur on the scale Xp- As a result, the quasi- 
particle description remains valid when the elastic mean free path, I ^ Xp. 
This translates to Ep ^ /i/r, where r is the elastic scattering time related 
to / by the Fermi- velocity. Multiple impurity scattering leads to the diffusive 
propagation of the quasi-particles for times » r or frequencies oj <^ 1/r, 
characterized by the quasi-particle diffusion constant, D. Fermi-liquid in- 
teractions between the quasi-particles, however, modifies D in a scale (fre- 
quency/temperature) dependent way, making the conductivity a, which is 
related to D through Einstein's relation a = {dn/dfj,)D, temperature de- 
pendent. Altshuler, Aronov and Lee showed [12, 13] that the corrections are 
localizing in the limit of weak quasi-particle interactions, so that the elec- 
tronic system becomes an Anderson insulator as the temperature T — >■ 0. 
Here, dn/d^ is the thermodynamic density of states which is not renormal- 
ized in the diffusion limit. (Experimentally, the compressibility is found to 
change sign in the strong-localization regime [14, 15], which is outside the 
scope of the diffusion theory.) 

In the discussion above the Fermi-liquid scattering amplitudes were con- 
sidered to be scale independent, as is the case in a clean Fermi- liquid. It 
was shown, however, that the disorder averaged amplitudes are in fact scale 
dependent at low energies, w < 1/r, due to the diffusive propagation of the 
quasi-particles [16]. For weak disorder (low resistance) the effective scatter- 
ing amplitudes were shown to become strongly correlated over long distances 
and eventually diverging at a very small but finite scale (temperature) T*. 
The enhanced correlation overcome localization so that the resistance p* at 
r* is finite. This interplay of diffusion and interaction controls the physics of 
dirty- Fermi-liquid down to the scale T*, below which the liquid is no longer 
describable by diffusing Fermi-liquid quasi-particles. 

In contrast to the weak disorder case, it was shown recently [3] that the 
situation is very different when the disorder (resistance) is increased. For not 
so weak disorder the Fermi-liquid correlations were shown to be suppressed 
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at long distance, so that the dirty-Fermi-liquid becomes an Anderson-like 
insulator at low temperatures. It was also established in Ref. [3], in the 
framework of an effective large- dirty-Fermi-liquid theory, that a QCP at 
a finite critical resistance separates the metallic and the insulating phases, 
implying the existence of a MIT in two dimensions. The theoretical details 
of these results and their experimental tests are briefly discussed below. 

2.1 Theory: Background 

The minimal model that describes the problem of quasi-particle diffusion in 
the field of impurities and interactions is the non-linear sigma model [16] 



S[Q] = 5dis + Sint 
•S'dis = 



^ I d'rTv[D{VQf-4z{eQ)] 



int 



+— / d'r Q(ri + r2 + rc)Q 



(la) 

(lb) 

(Ic) 



The functional S\Q] describes interacting quasi-particles with energies 
e < l/r. (e is the energy matrix whose elements involve the energy of the 
quasi-particles.) Note that the region of higher energies l/r < |e| < Ep, 
containing Fermi-liquid renormalizations of the clean liquid are incorporated 
as the input parameters into this theory. The Q-field here is an auxiliary 
matrix field describing the impurity scattering. It satisfies the constraints: 

= 1,Q = (5^, and Ti Q = 0. The fluctuations of Q in the particle- 
hole (diffusons) and particle-particle (cooperons) channels are described by 
the diffusive propagators 'D{q,ujn) ~ ^/{Dq^ + zcon) (we use Matsubara fre- 
quencies). These propagators capture the diffusive evolution of the charge 
and spin density fluctuations at large times and length scales, i.e., small fre- 
quency and wave- vectors. The parameter z is the frequency renormalization 
parameter [16, 17], which can be interpreted as the renormalization of the 
density of states of the diffusion modes (z = 1 for non-interacting electrons), 
and as a result controls the thermodynamic properties of the dirty-Fermi- 
liquid. The cooperons capture the effects of quantum interference, which at 
the perturbative level leads to the well known weak-localization corrections. 

The quasi-particle interactions are characterized by the Fermi-liquid am- 
plitude matrices, Fi and in the diffuson channel and, Fc, in the cooperon 
channel, respectively. The details of their matrix structure in the space of 
spin, frequency and momentum are not shown here, they can be found in 
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Ref. [1]. (See Ref. [18] for their valley dependence.) The magnitudes of the 
amplitudes, Fi and are related to the standard zeroth order harmonics 
of the Fermi-liquid amplitudes Fq''^ as Fi^2 = — -^o'"/(l + -^o'*^)- 

The scaling of the parameters, D, z and Fj's (i = 1,2), together describe 
the transport and thermodynamic properties of the dirty-Fermi-liquid. Spe- 
cial care is, however, required when Coulomb interactions are present [16]. 
The Fi amplitude in this case includes amplitudes of the kind which can be 
separated by cutting the statically screened long-ranged Coulomb line once. 
They are usually denoted as Fq. This distinction is important because the 
polarization operator which is irreducible to cutting a Coulomb line does 
not include Fq. Conservation laws then require that z = 2N{Tq + Fi) — F2, 
where is the number of valley degrees of freedom [16, 18]. 

As discussed in the last section, the diffusion constant and the Fermi- 
liquid scattering amplitudes become scale dependent in a dirty-Fermi-liquid. 
It has been shown that renormalization group (RG) theory applied to the 
dirty-Fermi-liquid model given in (1) is able to capture the scale dependence 
originating from the interplay of disorder and interactions to all orders in the 
interaction amplitude, making it the most promising analytical tool available 
to understand the physics of disordered systems. Pedagogical reviews of 
the RG theory can be found in Refs. [1, 2]. For more recent advances, see 
Ref. [19]. 

The scaling variables are described in terms of the dimensionless param- 
eters, the resistance p = l/2N{2'a'^uD), and the interaction F2 and z. It 
follows, however, from the structure of the action in (1) that the variables 
can be expressed as p and 72 = F2/2; which does not contain z explicitly. 
In terms of these variables, the equations for p and 72 form a closed set of 
equations independent of z. In general, they take the form: 

P/5p(p,72) (2a) 

P/572(P,72) (2b) 

P/3.(/0,72) (2c) 

The scale ^ = log(l/rr), where to logarithmic accuracy 1/r is used as the 
upper cut-off. Note that the diffusive contribution involving the parame- 
ter 71 is universal when the dynamically screened Coulomb interaction is 
involved, i.e., when 70 = Tq/z is included. Hence it does not appear as 



amp 
drf2 
dlnz 
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an explicit parameter in these equations [12]. It is a well established re- 
sult with great importance for the general structure of the theory that the 
static amplitude Fq does not acquire diffusion corrections and is therefore 
not renormalized. When combined with the conservation laws discussed 
above it follows that d{2NTi — T2) = dz, from which one may extract the 
renormalization of Fi [16]. The above discussions imply that the electronic 
transport (p) and thermodynamic (involving z) properties are describable 
by a two-parameter scaling theory involving p and 72 only. (For repulsive 
interactions, the amplitude 7c scales to a resistance dependent value for fi- 
nite N and hence is not considered here. The situation in the N 00 limit 
is different and will be discussed later.) 

The /3- functions are multiplied by a factor p to emphasize that the sought 
after corrections appear as a result of the underlying diffusive propaga- 
tors. The complete form of the /3-functions are unknown. The general 
approach, however, is to expand the functions in a power series in p as 
/3(/9, 72) ~ /3i(72) + ph{l2) + ■ ■ ■ ; such that for each power of p the full 
dependence on 72 is obtained. This is possible, in principle, because the 
maximum number of allowed interaction vertices are limited (before the 
ladder sums are included) by the number of momentum integrations involv- 
ing the diffusive propagators, and since each propagator gives a factor of 
1/Z) ~ /?, for a given order in p only a finite number of vertices contribute. 
So far, only the first order (one-loop) results have been obtained exactly [16]. 
Such a solution has the obvious limitation that it is applicable only when 
/9 <C 1 and is not applicable near /? ~ 1 where the MIT is experimentally 
found to be located [4]. Therefore to approach the MIT, the disorder has 
to be treated beyond one-loop (at least to two-loop order) while adequately 
retaining the effects of interaction. Exact results to order p^ have not been 
obtained so far. It was shown, however, that certain simplifications that 
come with standard large-A^ theories can be exploited to obtain the two- 
loop results [3]. The RG equations to two-loop order within the large- 
approximation showed, for the first time, the existence of a QCP that de- 
scribes the MIT separating the metallic from the insulating phase in 2D. 
Analysis of the QCP revealed a non-Fermi-liquid fixed point with diverging 
specific heat and spin-susceptibility. The one and two-loop RG equations 
are presented below and comparison with experiments are discussed. 

2.2 Finite-iV: One-loop 

The theoretical details of the evolution of p and 72 in this region with p ^ 1 
were discussed in detail in Ref. [5]. For completeness we give the equations 
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and only discuss its salient features here. The RG equations generalized to 
include equivalent valley degrees of freedom are: 

dlnp 



drf2 _ 
dC ~ 2 



p [N + l- (4iV2 - 1)^(72)] (3a) 
"(I + 72)' (3b) 



__ = ^[_i + (4iv2-l)72] (3c) 

where $(72) = (1 + I/72) ln(l + 72) — 1- The factor (4A^^ — 1) corresponds 
to the number of spin-valley triplet channels; the factor N in Eq. (3a) corre- 
sponds to the weak- localization corrections, which is enhanced due to the 
extra degrees of freedom. The factor of one in (3a) is the contribution of the 
long-ranged Coulomb singlet-amplitude, which after dynamical screening is 
universal; it should be emphasized that the factor of one arising in equation 
(3b) for 72 has the same origin, and as a result, setting the initial value of 
72 to zero does not switch the interactions off. 

The following salient features are to be noted: (1) while the amplitude 72 
increases monotonically as the temperature is reduced; (2) the resistance, as 
a result, has a characteristic non-monotonic form changing from insulating 
behavior [dp/dT < 0) at high temperatures to metallic behavior {dp/dT > 
0) at low temperatures. The change in slope occurs at a maximum value 
Pmax at a temperature T = Tmax, neither of which are universal; (3) the 
corresponding value of the amplitude 72 at Tmax is, however, universal at 
the one-loop order, depending only on A^; for = 1, it is 2.08, whereas for 
A^ = 2, it has the considerably lower value 0.45; (4) it follows from the form 
of the equations that the temperature dependence of p(T)/pmax and ^2{T) 
are universal when plotted as functions of Pmax ln(r/Tinax)- 

As explained in the introduction, the theory described here is applicable 
only in the diffusive regime. It is critical that the diffusive regime (temper- 
ature and density range) is identified carefully before the scaling is tested. 
In Ref. [6] a method that takes advantage of the B and T dependence of 
the magnetoconductance, a{B,T), in parallel field was exploited to identify 
the diffusive regime. It is known that while a{B,T) scales as (B/T)'^ in 
the diffusive regime [20, 21], in the ballistic regime it scales as B'^/T [22]. 
This distinction was used to reliably identify the diffusive regime in a MOS- 
FET sample. The experimental details and the results of the comparison 
with theory can be found in Ref. [6], where for the first time the scaling of 
the interaction amplitude was established. Remarkably, the parameter 72 at 
T = Tinax was found to correspond to 0.45 as predicted by theory for N = 2. 
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2.3 Two-loops and the MIT in the hmit iV ^ oo 

It is instructive to analyze Eq. (3) for ^ 1. In this case it is easily 
seen that the value of the interaction amplitude 72 at T = Tmax; i-e-, when 
dp/dT = in Eq. (3a), scales as -i2{Tra-^^; N) ^ 1/2N for large N. Further- 
more, because the presence of A'^- valleys increases the number of conducting 
channels, the resistance scales as /O ~ \/N. This can be explicitly taken 
into account by defining the maximum resistance, which is the only free 
parameter at the one-loop level, as Pmax = ^/N. Defining new variables 
t = p{T)/p^^s.y^ and 62 = 12{T) h2{T^a.y.] N) , and keeping in mind that 
» 1, Eq. (3) can be expanded in powers of 1/N: 

^ = t[l-e2 + 0{92/N)] (4a) 
dft 

^ = t[l + 0{e2/N)] (4b) 

^ = t[e2-0{l/N)] (4c) 

After rescaling the temperature scale becomes 77 = ln(rmax/r)- It is now 
possible to take the limit N ^ 00, keeping t and 02 finite. Note that the 
parameter t corresponds to the resistance per valley, t = 1/{2it)'^vD, with 
being the density of states of a single spin and valley species. In the 
infinite- A^ limit Eq. (4) reduces to 



dint 
drj 
dB2 
drj 
din 2; 
drj 



t{l - 62) (5a) 
t (5b) 
162 (5c) 



The form of Eq. (5) illustrates the fact that in the N ^ 00 limit the maxi- 
mum number of interaction amplitudes that can appear are limited by the 
power of t. This is to be compared with the function $(72) appearing in 
Eq. (3). Clearly, this is one of the main simplifications provided on taking 
the N ^ 00. 

Taking advantage of this simplification, the one-loop results above were 
extended to two-loops, i.e., order in Ref. [3]. It is easy to see that, unlike 
Eqs. (3), the solutions of Eqs. (5) are valid up to T = 0, at which point 
02 ^ 00 and t ^ 0, providing an asymptotically exact theory all the way 



9 



down to T = 0. (As discussed earlier, 72 for finite N diverges at a certain 
finite temperature T* befow which the RG equations are no fonger valid.) 
It is therefore left to analyze the large t behavior at two-loop order. 
The RG equations in the N ^ 00 limit are given as [3]: 



dint 
drj 
dB 
drj 
d\n.z 
dri 



t{a-Q) + t'^{l-a-ae + ctQ'^) (6a) 



a 



t{l + a + aQ) - At' (1 + a)G + -G^ + ceO'^ (6b) 

te (6c) 



where q = (5 — 7r^/3)/2 and ce = (1 — 7r^/12)/2. The parameter a takes 
values and 1, depending on if the cooperon channel is absent or present, 
respectively. In the large- limit because rescattering in the cooper channel 
is not accompanied by factors of N such processes are irrelevant and hence 
6c = N^c appears on equal footing as 62- (Violation of time-reversal by a 
magnetic field for example will, however, suppress the contribution of 7c.) 
The amplitudes 62 and 6c always appear in the combination Q = 62 + a6c- 

It is clear from Eq. (6b) that for large enough resistance dQ/drj < 0. 
Suppression of the 6 amplitude leads to allows t to grow to large values 
(insulator) as the temperature is lowered. This is a clear departure from the 
metallic behavior observed for small t. The insulating and metallic phases 
are then separated by a fixed point given at {tc, @c) where dlnt/dr] = and 
dQ/drj = 0. The quantum critical nature of the fixed point is apparent from 
the equation for z which at {tc, &c) diverges as z ~ l/r*'^®'^. The parameter 
z, which corresponds to the frequency renormalization, can be interpreted 
as the density of states of the underlying diffusion modes. Therefore it con- 
trols the behavior of the specific heat ~ {zu)T. Hence, divergence in z 
implies a divergence in the specific heat constant C^/T signaling the break 
down of the Fermi-liquid behavior at the critical point. Similar divergence 
is expected in the Pauli spin susceptibility also. Since the interaction pa- 
rameter G is finite at the critical point, the divergence in the Pauli spin 
susceptibility is not related to any Stoner-like magnetic instability. 

The feature that can be compared with experiments comes by observing 
in Eqs. (6) that unlike in the metallic phase where dt/dij < and dQ/dij > 0, 
in the insulating phase the opposite behavior is obtained, namely, dt/dr] > 
and dQ/dr] < 0. In Ref. [6], precisely such a scaling was confirmed experi- 
mentally. 
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3 Conclusions 



We believe thereby that the two-parameter theory of the dirty- Fermi-Uquid 
reviewed here captures both quantitatively (for small disorder) and qualita- 
tively (for large disorder) the physics of the MIT in 2D. 
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